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Particle production due to a quantized, massless, minimally coupled scalar field in two-dimensional 
flat spacetime with an accelerating mirror is investigated, with a focus on the time dependence of 
the process. We analyze first the classes of trajectories previously investigated by Carlitz and Willey 
and by Walker and Davies. We then analyze four new classes of trajectories, all of which can be 
expressed analytically and for which several ancillary properties can be derived analytically. The 
time dependence is investigated through the use of wave packets for the modes of the quantized field 
that are in the out vacuum state. It is shown for most of the trajectories studied that good time 
resolution of the particle production process can be obtained. 

PACS numbers: 03.70.+k, 04.62. +v, 04.60.-m 



I. INTRODUCTION 



As the simplest theoretical manifestation of the dynamical Casimir effect (DCE), the moving mirror model of 
DeWitt p], and Davies and Fulling [3] describes the disturbance of a field by an accelerated boundary, which 
results in both particle production and a flux of energy. As the mirror model matured [3HH] > it became apparent that 
accelerating boundaries could be used to understand entropy production [§1 [TU] , the relationship between particles and 
energy [TT], and thermodynamical paradoxes [HHH]- The DCE has the potential to be measured [13 US]- Indeed, 
one experiment claims to have effectively measured the DCE using a superconducting quantum interference device 
(SQUID) that acts as a moving mirror [T7]. In another recent set of experiments [TH] an analogue of the DCE has 
been observed in the case of a Bose-Einstein condensate. 

One aspect of the moving mirror model that has been largely unexplored is the study of the time dependence of the 
particle production process. The time dependence of the stress-energy tensor for the quantized field has been worked 
out. However, the stress-energy contains vacuum polarization effects along with particle production effects, and in 
most cases there is no clear way of separating the two. In a flat-space background the Bogolubov transformation 
between the in and out vacuum states can be used to accurately describe the particle production process. However, 
the resulting particle frequency spectrum retains no information regarding the time dependence of creation. 

One way to explore the time dependence of particle production is through the use of wave packets. Hawking |19j 
made use of such packets to describe the late-time behavior of black hole radiation when a black hole forms from 
collapse. In the context of a moving mirror, wave packets have been used by Dorca and Verdaguer [20] for a specific 
class of trajectories that generate a thermal spectrum at late times. 

In this paper we explore the use of wave packets as a means of obtaining information about the time dependence 
of particle production due to accelerating mirrors in (1+1) dimensions. We do so using two previously investigated 
trajectories and four new ones. In each case we compare the particle creation results with the stress-energy of the 
quantum field. We work with a massless minimally coupled free scalar field. The packets are obtained by integrating 
the modes of the quantum field over specific frequency ranges using a parameterized weighting function that creates 
a packet whose amplitude is largest near a particular time that is related to the value of one of the parameters. By 
computing the Bogolubov transformation using packets for the modes in the out vacuum state it is possible to obtain 
an expression for the number of particles produced in various frequency and time intervals as a function of time [21) . 
There is a fundamental uncertainty principle involved in working with the wave packets such that small frequency 
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bins lead to good resolution in frequency and poor resolution in time, and vice versa. We explore the effects of this 
uncertainty relation for the trajectories chosen. 

To facilitate the investigation we have restricted our attention to mirror trajectories for which the Bogolubov 
components can be computed analytically. Only a small number of classes of trajectories have previously been 
considered for which it is possible to analytically compute the Bogolubov coefficients. These include the trajectories 
studied by Carlitz and Willey [7], which are designed so that there is constant flux of energy, and the class of trajectories 
studied by Walker and Davies |22| . which involve a mirror that starts at rest in the infinite past, accelerates, and ends 
at rest in the infinite future. In this paper we introduce four new trajectories for which the Bogolubov coefficients 
can be computed analytically. Three of these involve mirrors that start at rest in the infinite past. In one case the 
mirror trajectory is also asymptotically static in the future, in another the mirror's speed is asymptotically constant, 
and in the third the mirror's speed approaches the speed of light. In this third case the trajectory nevertheless ends 
up at future timelike infinity instead of being asymptotic to a null trajectory. A fourth class of trajectories begins in 
the same way as those studied by Carlitz and Willey but then, instead of becoming asymptotic to a null ray at late 
times, approaches a constant velocity and thus becomes inertial. 

In the only previous use of wave packets for the moving mirror model that we are aware of, Dorca and Verdaguer |20j 
studied the asymptotic form of the trajectory originally discussed by Davies and Fulling [HE]. Their procedure involved 
wave packets for modes in both the in and out states. This technique allowed them to obtain a finite spectrum even 
though the total number of particles produced by the mirror in that model is infinite. 

Here we use wave packets only on the modes that approach future null infinity to the right, J^^ , and are in the out 
vacuum state. The idea is to model what a particle detector at a large distance from the mirror would see. It turns 
out that one can compute either the packets first and then the Bogolubov transformation to determine the number 
of particles in a given packet or one can do the computation in reverse order. We do the latter and first compute 
the exact Bogolubov transformation for the modes and integrate over frequencies using the appropriate weighting 
function [5T] to obtain the Bogolubov coefficient for a packet. To obtain the total number of particles in the frequency 
range for that packet we then integrate over the entire frequency range u/ of the modes in the in vacuum state. For 
two of the trajectories considered there is an infrared divergence in the number of particles created, which manifests 
as a divergence in the number of particles in those packets that include modes of arbitrarily small frequencies. Any 
real detector would have an infrared cutoff, so we eliminate this divergence through the use of a low frequency cutoff 
in the computation of the wave packets. That is, for these trajectories, we do not consider packets that include modes 
in the out state with arbitrarily small frequencies. In two other cases the Bogolubov transformation results in an 
infrared divergence in the total number of particles created but there is no corresponding divergence in the number 
of particles in the packets with small frequency modes. In the other two cases there are no infrared divergences and 
the total number of particles is finite. 

The outline of this paper is as follows. In Section [TTJ we review the relevant aspects of the moving mirror model 
in (1+1) dimensions in some detail. In Section III we discuss both the Carlitz- Willey [7] and Walker-Davies [22] 
trajectories. For the Carlitz- Willey trajectory we first review some previous results and provide an analytic expression 
for the trajectory. We then use wave packets to show that the spectrum is time independent with a Planck character. 
We analytically compute the spectrum at a given time for wave packets of various frequency widths to investigate the 
distortion of the spectrum. For the Walker-Davies trajectory we review some of the analytic results found previously. 
Unlike the Carlitz- Willey case, it does not appear possible to calculate analytically the expected number of particles in 
the wave packet modes. In Section [TV] we discuss four new trajectories and for each compute the relevant Bogolubov 
coefficients, the stress-energy tensor, the number of particles produced at a given frequency, and for one trajectory, 
the total number of particles produced. Section [V] includes a comparison of the time dependence of the stress-energy 
tensor for the quantum field and the number of particles produced, where the latter is investigated through the use 
of wave packets. The difficulties encountered in obtaining simultaneously good time and frequency resolution for the 
number of particles produced are discussed, as well as the divergences that occur for some trajectories in the number 
of particles produced and the energy of those particles. Our conclusions are given in Section VI Throughout this 
paper units are used such that H = c = 1. 



II. BACKGROUND 



The moving mirror model in (1+1) dimensions consists of a quantized massless scalar field in flat space that obeys 
Dirichlct boundary conditions on a perfectly reflecting (mirror) boundary. The scalar field satisfies the wave equation 



□$ = 



(2.1) 
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In this paper we always expand the held in terms of mode functions that are parameterized by the frequency u>. 
Denoting them for the moment as <j> u , they obey the equation 



{-d 2 t +d 2 x ) </> u = -a u d v < l 



o, 



with 



u = t — x , 

V = t + X . 

The general solution is 

0a, = g(v) + h(u) , 

with g and h being arbitrary functions. The mode functions are normalized using the scalar product 



■l I dZn^ 



<p! (x) df, 4>* 2 {x) 



(2.2) 



(2.3a) 
(2.3b) 



(2.4) 



(2.5) 



Here E is any Cauchy surface for the spacetime, n M is a future-directed unit normal to that surface |21j , and we adopt 
the usual notation for the derivative acting to the right first and then to the left with a minus sign 23J. For this 
scalar product the canonical relations hold, 



{(j) u) {x),<j) ul ,{x)) = -(4>*(x),<f>Z,(x)) = S(ui - u') , 



(2.6a) 
(2.6b) 



For Minkowski space with no boundaries we can choose the normalized modes 

1 



\/47ru; 
1 
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Then 



dui 
o V '4ttu 



(2.7a) 
(2.7b) 

(2.8) 



with a uu , a uv , a\j U , and aj^ being the usual annihilation and creation operators. 

If there is a mirror with trajectory z(t), so that at any time t the mirror is at the position x = z(t), then the 
spacetime effectively has a boundary. In this paper we only consider solutions to the mode equation that are to the 
right of the mirror and incorporate reflection from the mirror's surface. We also only consider mirror trajectories that 
begin at past timclikc infinity, i~ . In this case past null infinity, £~ , only consists of the surface u — — oo. This 
is a Cauchy surface. If the mirror trajectory ends at future timclikc infinity, i + , then future null infinity, J^ + , only 
consists of the surface v = oo and this is a Cauchy surface. But if the trajectory is asymptotic to the null ray v — Vq, 
then has two parts, l/£ and ^ r £ l ", using the notation of Ref. [7]. The surface lies at v = oo and J?^ consists 
of the part of the surface u — oo which goes from v = Vq to v = oo. Taken together they also provide a Cauchy 
surface. 



It is useful to evaluate the scalar product (2.5 1 using the Cauchy surfaces for J? and y + . It can be shown (see 
e.g. [24]) that for J~ 



and for 



'1,02 J = 



<pi(u — —oo,v) d v (f>2 {u — —oo,v) 



dv , 



(2.9) 



4>i (it, v = oo) d u (f>2 (v., v = oo) 
0i (u = oo, v) d v 02 ( u = oo j v) 



du 
dv . 



(2.10) 
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If Dirichlet boundary conditions are imposed on the scalar field, then the mode functions 4>u must vanish at the 
location of the mirror. To quantify this, it is useful to introduce functions u m (t) and v m (t) that give the values of u 
and v at the location of the mirror at a given time t. Thus 



u = u m (t) = t — z(t) 
v = v m {t) = t+ z(t) 



(2.11a) 
(2.11b) 



We can invert the first equation to get t via a function t m (u) or we can invert the second to get t via a function t m (v). 

In the presence of a mirror we can consider cither the mode functions that arc positive frequency at , and thus 
correspond to the in vacuum state, or the mode functions that are positive frequency at and thus correspond to 
the out vacuum state. The modes that are positive frequency at are 



1 



\/4 



7TU/ 



e -iu'v _ e -iu'p(u) 



(2.12) 



Substitution into Eq. (2.6) using Eq. (2.9) shows that the normalization is correct. For these mode functions to vanish 
at the mirror we must have v — p(u) at the location of the mirror. If we invert Eq. (2.11a) above to find t = t m (u) 



and then use the definition (2.11b), we find that 



p(u) = t m (u) + z(t m (u)), 



(2.13) 



which fixes the function p(u) . 

In a general left-right construction (for mirrors that have a horizon at vq), there are two sets of mode functions 
that are positive frequency at y + . One set, which we will denote as (f>^, are nonzero at .JP^ and zero at J^j^] The 
other set, </>^, are zero at J't and nonzero at J^. The former are given by 



iR, out 



1 



V < Vq 



(2.14) 



Mirrors that are asymptotically incrtial in the future have vq = oo. Note that all trajectories that begin at past 
timclike infinity, i~ , (the only type we consider here) span the range — oo < u < oo. Substitution into Eq. ( |2.6[ ) using 
Eq. (2.10) shows that these modes are normalized correctly also. Again, for these modes to vanish at the mirror, 
we must have u = f(v) at the location of the mirror which is the inverse relation to v — p(u) and an equivalent 
representation of the mirror trajectory in u, v coordinates. If Eq. (2.11b) is inverted to find t = i m (v) and the 
is used, then one finds that 



definition (2.11a 



f(v) = t m (v) - z(t m (v)) 



(2.15) 



which fixes the function f(v). 

There are no other modes if vq = oo. But if the mirror's trajectory is asymptotic to the null surface v — Vq, then one 
must also include the set of modes <fi^ that reach J^^~ and never interact with the mirror. Substituting into Eq. (2.6) 



and using Eq. (2.10) one finds that 



(2.16) 



To further examine the behavior of these modes it is useful to work with a specific trajectory. This has been done by 
Carlitz and Willey [7] for a trajectory with a future horizon. We study other aspects of this trajectory in Sec. 



Ill 



The usual procedure for calculating interesting observable quantities, such as the energy or particle number, starts 
with the choice of an appropriate trajectory, z(t). Then either the function p(u) or its inverse f(v) is found using the 
procedures described above. A key aspect of these procedures involves the solution of the relevant, and sometimes 
transcendental, function inversions. This requirement has made it difficult to find trajectories that allow both a fully 
analytic description of the mirror's motion and an analytic calculation of the associated Bogolubov coefficients (see 
below). 

The function p(u), commonly called the ray-tracing function [26], characterizes the mirror trajectory and is incor- 
porated in the modes, the two-point function, the energy flux, and the correlation functions. The trajectories and, 
where known, the ray-tracing functions for the mirrors considered in this paper are given in Table [I] 



1 The left and right coefficient formulation lies at the crux of the calculations in Ref. [J| and despite the call for more attention [25], the 
construction has been under-utilized. 
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TABLE I. Some classes of trajectories, z(t), and the ray-tracing functions p(u) and f(v) associated with them. 



Trajectory p(u) f(v) 



Static 


z — 





P = 


u 


/ = 


V 


Constant velocity 


z — 


—Vot 


P = 


1 -"Q ? , 


/ = 


l-VQ 


Uniform acceleration 


z — 


K -1 - Vk- 2 + t 2 


P = 


U 

l-\-K,U 


/ = 


V 

1 — K.V 


Carlitz-Willey 


z — 


-t-±W{e~ 2Kt ) 


P = 


_± e -™ 

K 


/ = 


-ilog(-K«) 


Walker-Davies 


t = 


-z±AVe- 2z / B - 1 










Arctx 


z — 


-itan-^e'") 










Darcx 


z — 


-isinh-^e" 4 ) 










Proex 


z — 


-\w{e^) 


P = 


u- \ W{2e pu ) 


/ = 


v+ 2 p e^ 


Modified Carlitz-Willey 


z — 


-W[e-^/(l + a)]/ K 


P 


c7ii — \e~ KU 


/ = 


I + \W{e~ VKia jo 



A. Bogolubov transformations 



One way to describe the particle production that arises in the presence of an accelerating mirror is to use the 
Bogolubov transformation. The positive frequency modes at </>™, form a complete set and one can expand 

modes at ^ + in terms of them, 



with J representing either R or L. Using the relations (2.6 1 one finds 



^-(^,C'), 



/?1' = -(<^,C<*). 

The field <& expressed in terms of the mode functions can be represented in either of two ways, 



duJ 



m / m , 



tot. 



/"DO 

7 JO 



Using 6^ = one finds 



6 J ,= 



dJ 



./ s* inf 



(2.17) 



(2.18a) 
(2.18b) 



(2.19a) 
(2.19b) 

(2.20) 



If the field is in the in vacuum state specified by the positive frequency modes at ^ , we can use the operator 
= (b'^Ybfj to compute the average number of particles with frequency u that reach <#j~ , 



{N J U ) = {Q in \N J M = / dJ\f3L,\ 2 



The expectation value of the total number of particles that reach J?j~ is 



/' oc poo 

(N J ) = (Oml^lOta) = / dw / du/ |/3, 
Jo Jo 



>\qJ |2 



(2.21) 



(2.22) 



Since we are primarily concerned with the number of particles that reach J'^ , we will focus on the computation of 



P% u ,. If the Cauchy surface J~ is used, then Eq. ( |2.18b[ ) along with Eqs. (2.9| >, ( |2.12[ ), and ( |2.14[ ) gives 

1 



(2.23) 
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If the Cauchy surface is used, then one similarly finds 



B R , = 



du e-*"»-*» | u 



,dp(u) 
du 



(2.24) 



These expressions are of cour se eq uivalent. 

It is possible to write Eq. (2.24) in terms of a time integral over a function of the trajectory z(t). Since p(u) is a 



fixed function of u and u ranges from 
to change variables, one finds 



-oo to +oo, one can substitute u m (t) for u in Eq. (2.24). Then using Eq. (2.11a 



oR 



dt e -*"+*+-*"-*(*) (u+z{t) - U-) 



(2.25) 



where oj + 



and 



Note that if we consider an inertial trajectory it is easy to show that f3 R , = 



Thus, as expected no particles are produced when the mirror does not accelerate. 

It is not hard to show that the Bogolubov coefficient a uu i may be obtained from the above expressions for fj uu i by 
letting {J — > — uj' everywhere in the expressions for B^', except for the factor \/\f luuj' which must remain unchanged. 
Finally, we note that if the trajectory is initially inertial and the acceleration does not continue forever, then the total 
energy produced is finite and given by the following sum over the quantum modes [llj 



E„ 



oj(N u )du! 



(2.26) 



B. Stress-energy tensor 



The renormalized stress-energy tensor for the massless, minimally coupled scalar field was computed in terms of the 
function p(u) and in terms of the trajectory z(t) by Davies and Fulling [2J. They found that the energy flux produced 
by the mirror as a function of u is given bj|^] 



1 

247 



P_ 
p' 



(2.27) 



where primes indicate derivatives with respect to u. Their expression for the energy flux in terms of z(t m (uj) is 

'i'(i 2 - 1) - 3zz 2 



I2n(z- l) 4 (i + l) 2 



(2.28) 



t=t m (u) 



where the dots refer to derivatives with respect to t. Eq. (2.28) is equivalent to Eq. (2.27) evaluated at the surface of 
the mirror. In either case, it is easy to show that for an inertial trajectory (T uu ) — 0, as would be expected when the 
scalar field is in the vacuum state. 

It is also possible to write Eq. (2.28) in terms of the time derivative of the proper acceleration^] 



The result is 



a = z/{l-z 2 f' 2 . 
a \/l + z 



12tt (1 - z)3/ 2 



(2.29) 



(2.30) 



The overall negative sign implies that on the right hand side of the mirror a flux of negative energy is given off if 
the change in acceleration of the mirror is towards the right and a flux of positive energy is given off if the change in 
acceleration is towards the left. 

For trajectories that are asymptotically inertial in the limits t — > ±oo, p(u) — > c\u + cq for some constants c± and 
cq. For trajectories considered in this paper, a finite amount of energy reaches (except for the Carlitz-Willey or 



2 Some other components are (T M > = (T tx ) = (T xx ) = (T uu ). 

3 The proper acceleration is the acceleration in the instantaneous rest frame of the mirror. Note that the time derivative is in the inertial 
frame, not the rest frame of the mirror. 



7 



modified Carlitz-Willey classes). Because the flux (2.271 is only a function of u and therefore does not fall off at J^, 
the total amount of energy E st that reaches can be obtained by integrating over u, 



(T uu ) du 



(2.31) 



Walker has a proof that E st — E qs provided that the mirror is asymptotically inertial in both the past and future, 
i. e. cc(±oo) = 0, and the velocity towards J?^ never reaches the speed of light, i(±oo) ^ 1. Substituting Eq. (2.27) 
into (2.31) and integrating by parts one finds 



E, 



1 

48^ 



oo / // \ 2 

p 



p' 



du 



(2.32) 



so long as the surface terms with p" jp' vanish asu-> ±oo. The result can also be written in terms of a time integral 
over a function of t he tra jectory of the mirror by letting u — > u m (t), using Eq. (2.11a) to invert to find t = t m (u), 
and then using Eq. (2.28). The result is 



E^. — —-. 



1 

127T 



a\l - r dt 

oo 



Integrating by parts gives another expression 

1 



E s t — 



12tt 



a z (l + z)dt 



1 

127T 



(l + i) 2 (l-i) 3 



dt 



(2.33) 



(2.34) 



Notice that this last form masks the dependence on the time derivative of the proper acceleration. 



C. Wave packets 

Another way to investigate particle production is to use wave packets [IHJ [H] ■ An advantage of this approach, as 
discussed in the Introduction, is that one can study the time-dependent aspects of particle production. 

A wave packet, <f»j n , can be constructed from <j) u by integrating over a finite range of frequencies with a particular 
weighting function so that [2"T] 



1 



du e 2 ™ n/e <j> u . (2.35) 



Here n takes on integer values and j takes on nonnegative integer values. Substituting Eq. (2.14) into Eq. (2.35) 



and noting that the first term does not contribute due to rapid oscillations in the limit v — > oo, one can see that the 



integral is largest for values close to u = 2irn/t. It is clear from (2.35) that the value of j is related to the frequency 
of the modes in the packet with (j + l/2)e giving the frequency at the center of the range and e giving the width of 
the range. When the weighting functions are applied to the modes e~ tulu / VAttoj , the resulting wave packets form a 
complete and orthonormal set. 

One can use the scalar product to construct the Bogolubov coefficients that correspond to the wave packets |21j . 
As mentioned previously, we are concerned in this paper with the particles that reach ■ In that case 

/9fn, w ' = -(0f„,C'*)- (2-36) 

It is possible to obtain these wave packet coefficients directly from the coefficients f}^ , by using the same weighting, 
integrating over frequency, and swapping the order of integration 



1 R = — 



duj e 2 ~ n / e ^, . (2.37) 



The average number of particles produced for given values of n and j is 

/>oo 

(JV«)= / <(u/|<3»,„,| 2 



This quantity gives the average number of particles that reach J 2 ^ in the frequency range je < u) < (j + l)e and 
in the approximate time range (27m — ir)/e < u < (27m + 7r)/e. It can be used to estimate the average number of 
particles that a detector would see in this frequency range if it was turned on during the above time period near an 
event centered at some x and some large v. Thus computation of (Nj n ) for a range of values of j and n allows one to 
construct the evolution of the spectrum of the produced particles in time, to the extent allowed by the uncertainty 
relation, as it would be seen by a scries of particle detectors spread out over a line of constant but large v. 

Note that one can also estimate the total energy of the particles produced by multiplying the number of particles 
in a given bin by the frequency at the center of that bin, 




(2.39) 



This estimate of the energy can be compared to the energy of particles produced E qs or the total stress-energy flux 
E st to test the accuracy of the wave packet description of particle production. 



III. TWO PREVIOUSLY-STUDIED MIRROR TRAJECTORIES 



In this section we examine two previously-studied types of mirror trajectories. One of these, developed by Carlitz 
and Willey [7, consists of a trajectory that has a future horizon at v — Vq — but no past horizon. The trajectory 
is designed to yield a constant stress-energy flux. The functional form of the trajectory allows many quantities of 
interest to be computed analytically. The second type is the class of trajectories studied by Walker and Davies [22]. 
For these trajectories, the mirror begins and ends asymptotically at rest. Thus the total number of particles produced 
is finite. In this case too, a number of quantities can be obtained analytically. For both of these types of trajectories 
we have extended the analysis by using wave packets to compute the spectrum of created particles. 



A. Carlitz- Willey trajectory 

In their paper [7J Carlitz and Willey point out that if the motion of the mirror is specified (in u and v = p(u) 
coordinates) by taking the ray-tracing function to be 



p(u) = e 

K 



(3.1) 



then a constant energy flux results. Substitution into Eq. (2.27) gives the energy flux in terms of the free parameter 

K, 



48tt 



(3.2) 



An implicit functional form of the trajectory in t and x = z(t) coordinates can be obtained by substituting Eq. (3.1 ) 
into Eq. (|2.13|), 



1 



t + z(t) = --( 



_ p — K,t~\-Kz(t) 



(3.3) 



Carlitz and Willey did not provide the explicit functional form for z{t). However, we find that it can be given as 

f 



z{t) = -t--W(e~ AKt ) 



(3.4) 



which involves the Lambert W function (also known as the product logarithm). A plot of this trajectory is given 
in Ref. [7J and is shown also in our Fig. [T] It is not difficult to show that i — > ±1 in the limits t — > =Foo and that 
z < for all time. The mirror trajectory begins at past timelikc infinity, i~ , and at late times approaches v = 0. 
Substitution into Eq. (2.29) gives the proper acceleration, 



a{t) 



(3.5) 



which is not constant, even though the energy flux is. 
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In Ref. [7] analytic expressions were found for the Bogolubov coefficients. In particular it was found thatQ 



oR 



2lo 



Thus 



3 R , 

UJUJ' 



-7TCJ/2K. 



1 



— iuj/h 



K 



1 



2ttku}' e 2 ™/ K - 1 



(3.6) 



(3.7) 



resulting in both infrared and ultraviolet divergences for the quantity (N^) in Eq. (2.21 ). This is not surprising given 
the fact that the acceleration of the mirror, while zero in the limit t — > — oo, is nonzero at any finite time in the past. 
The radiation is produced with a thermal spectrum and propagates to 



It is possible to compute analytically the expectation value of the wave packet number (Nj n ). First substitute 



Eq. (3.6) into Eq. (2.38), interchange the order of integration, and make the variable transformation y = lnu/. Then 
integrating over y gives a result that is proportional to <5(wi — 002)- Next integrating over UJ2 and setting uj\ = uj gives 



(Njn) 



1 



C7+1)* 



du 



1 



In 



1 



1 



- 1 



(3.8) 



There is a divergence in the particle count in the lowest frequency bin, j = 0. This result is similar to the infrared 
divergence found previously in (N^) and the reason for it is the same. Physically of course, one cannot measure 
particles of infinite wavelength. Thus in a real particle detector the lowest frequency bin would have a lower limit 
cutoff rather than extending all the way to w = 0. 

It is also evident that for this mirror there is no dependence in (Nj n ) on the parameter n and accordingly the 
average spectrum of particles recorded by a detector would be independent of time. This is almost certainly related 
to the fact that for this trajectory the energy flux is constant. 

To find effects of various frequency ranges e on the spectrum and to recover the Planck form for the spectrum in 
the limit that e — > 0, one can first write the expression in (3.8) in terms of ojj = (j + l/2)e and then expand it in 
powers of e with the quantity ojj fixed. The result is 



(Njn) 



1 



1 



6« 2 (e 2 ™i/" - 1)2 v ' 



(3.9) 



With e fixed, the second term approaches e 2 /(12w 2 ) for small ojj and behaves like 7t 2 £ 2 /(6k 2 ) for large values of ujj. 
Thus, for a given frequency width the deviation of the spectrum from the Planck form becomes more pronounced for 
smaller values of the central frequency Lij . 



B. Walker-Davies trajectory 

The Walker-Davies trajectory [22] is given by the relation 

t = -z±A\Je- 2z l B - 1 , (3.10) 

where A > B and where the plus sign is adopted for t > and the minus sign for t < 0. The trajectory is plotted 
in Fig. [T] Note that the curve is C°° in spite of the change in sign across branches for positive and negative values 
of t. The mirror begins at i~ at rest with t,z = —00. It accelerates to the right and then decelerates back to rest at 
t = z = 0. Then it first accelerates and then decelerates to the left, ending at rest at i + at t — 00 and z — — 00. The 
modes to the right of the mirror always end at ■ Therefore we drop the subscript R in what follows. 
Walker and Davies [22] calculated the stress tensor and found that r] 

_ B (x 5 + - 2A 2 X 3 ~ WA 2 X 2 - 3A 4 X - ^B) 

{uu) ~6ir {x 2 + 2B X + A 2 Y ' [ j 



4 Here we have adapted the expression given in [7] to the conventions we are using. 

5 This thermal spectrum is of course that of a one-dimensional black body. A similar physical manifestation of one-dimensional thermal 
radiation is the resistor. Discovered by Nyquist 1271 in 1928, a resistor in a lossless transmission line of great length in equilibrium at 
temperature T has thermal electric noise that is an analog of a black body in one dimension. 

6 Note the misprint in their expression for the numerator. 
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FIG. 1. Six analytically-known mirror trajectories. In the left panel the mirror trajectories are plotted in t and x coordinates. 
The right panel depicts the same trajectories in a Penrose diagram. The Carlitz-Willey trajectory, with k = 1, and the modified 
Carlitz-Willey trajectory, with k = 1 and a — 1/3, are shown as solid curves. The Carlitz-Willey trajectory gives rise to a 
horizon at v = 0, while the modified Carlitz-Willey trajectory ends at future timelike infinity The short-dashed curve 

indicates the Arctx mirror with y, = 1, which is static in the distant past and future. The long-dashed curve is the Darcx 
trajectory with v = 1 and asymptotic future velocity £ = —1/2. The Walker-Davies trajectory for A = 2, B — 1 is shown as the 
dot-dashed curve. The Proex trajectory with p — 1 is denoted by the dotted curve. We consider only the region of spacetime 
to the right of a mirror. Thus past null infinity on the left (not labeled) plays no role in our analysis. Similarly, there is an 
abbreviated portion of future null infinity on the left (^l ), but only in the case of the Carlitz-Willey trajectory. All of the 
other mirror trajectories begin and end at i~ and i + , respectively. 



where x is a parameter related to the null coordinate u by 



u 



= B\n(x 2 /A 2 + l)+ X . (3.12) 



By integrating the flux over all u, they found the total energy to be 



B 2 

Est = 48(A2 - 52)3/2 ' ( 3 - 13 ) 



They also were able to calculate the Bogolubov coefficients and found that 



l/WI 2 = ^ (-^— ) Binh(7ra;fl)|Jif,(r)| a , (3.14) 



7T 2 \Lj' + U> 

where q = — | + iuiB, r = A(uj' + ui), an d K q { r) is a modified Bessel function of the second kind. We have not found 
an analytic expression for (N u ) in Eq. ( 2.21| and therefore also do not have one for E qs . Furthermore, we do not 
have an analytic expression for (Nj n ) in this case and instead have computed it numerically. We defer until Sec. [V] 
discussion of these results. 



IV. NEW TRAJECTORIES 



In this section, we consider four new types of trajectories for which it is possible to analytically compute the Bogol- 
ubov coefficients /3 uu i. Having this analytic result in turn makes it feasible to compute numerically the expectation 
value of the particle number (Nj n ) associated with the wave packets. Trajectories of each type are plotted in Fig. W 
for specific parameter choices. The functional form of the trajectories and their corresponding ray-tracing functions, 
where known, are summarized in Table [Tj 
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A. Arctx mirror trajectory 



As seen in the previous section the Walker and Davies [22] class is composed of trajectories that begin and end 
at rest in the limits t — > ±00. They produce a finite number of particles and a finite amount of energy. Another 
trajectory with these features can be devised by taking 



z(t) 



— tan -1 (e"*) 
A 4 



(4.1) 



Here /i is a positive constant. Such a mirror starts at rest at x — in the infinite past and ends at rest at x = — 7r/2/^ 
in the infinite future. We refer to this trajectory as Arctx, drawn from Arctangent Exponential. 

The functional form of this trajectory is simple enough that a number of properties can be derived analytically. 
For example the proper acceleration is given by 



a(t) = j-l 



4 sech(/ii) tanh(/ii) 
[4 - sech(/xi)] 



3/2 



(4.2) 



Starting from zero in the t — > —00 limit, the proper acceleration is negative for t < 0. It first increases and then 
decreases in magnitude, before reaching zero at t — 0. For positive values of t it is positive, and again first increases 
and then decreases to zero in the limit t — > 00. A second important quantity, the stress-energy, can also be obtained. 
Substitution of p~T| into (|2.28[) gives 



yU 2 cosh[/it m (w)](— 5 — 2cosh[2^t m (u)] + cosh[4/ui m (w)]) 
3tt(1 - 2cosh[^ m (u)]) 2 (l + 2cosh[^t m (u)]) 4 



(4.3) 



The finite total energy would follow from integrating this flux over u, but it is simpler to substitute (4.1 1 into (2.34) 
to find 



(4.4) 



Next, it proves possible to calculate analytically the Bogolubov coefficients. Substituting (4.1) into Eq. (2.25) we find 

/W = 30 [giT(-m)T(-q) - g 2 T(m)T(q)] , (4.5) 

where 

(4.6a) 
(4.6b) 

: sin(7rm) 



/' 
1 

2^L 



(w' - w) 



,9o 



le 2 



yj 4m 2 



2ir/j, 

2-Plfl(l - ™, 1 



sin[7r(m + q)] 

91 = 2 ^ lM l -m,l-q,l-m-q, -l)e'^ m +^ 

92 = 2F 1R (l +TO, 1 + q, 1+TO + q, -1) , 

and the 2F1R are regularized hypergeometric functions. 



(4.6c) 

(4.6d) 
(4.6e) 



B. Darcx mirror trajectory 



Another analytically simple set of trajectories is found by setting 

z(t) = -^sinbT 1 (e 1 '*) , (4.7) 
v 

where v and £ are constants. To maintain future asymptotic inertial behavior, it is necessary that < |£| < 1. In this 
case, the mirror begins at rest and in the limit t — > 00 its velocity approaches — £. We refer to this set of trajectories as 
Darcx, short for Drifting Arc- Hyperbolic Sin of an Exponential. A specific example is plotted in Fig. [l] Remarkably, 
the Bogolubov coefficients (3 uu i and other relevant quantities for these trajectories can also be calculated analytically. 
Since the process of deriving them is identical to that discussed previously, we simply collect the results in Table [Tl] 
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TABLE II. Bogolubov coefficients and other useful information for the Darcx trajectories. 



Bogolubov coefficient /W = — )== -2 iu + 4 %^ r ( -^ +) r(ia + ) 

|2 £ 2 u/aj csch(7TOJ + )csch(7ra + ) 



with b+ = bco -\- aid , a+ = auj + bcu' , cu+ = - (w + a/ ) 
=i(l + 0,and6=i(l-0. 



Proper acceleration aft) = —7 ^ tttt- 

V V ' [l_( S 2_ 1)e 2„t].V2 

Rnertrv flnv IT \ - v - €xe" tm( " ) [2(g 2 -l) e 2 -^("> + i] 

Energy mix - i2¥ ( x+ee ^m(«))2[ K 2_ 1)e 2w m ( u )_ 1 



Total energy E st = ^ In ± 



l+i _ 3+g(3+2g) \ 
' « £(1+0 r 



C. Proex mirror trajectory 



Another interesting trajectory that is asymptotically inertial in the past, t — > —00, can be defined using the Lambert 
W function, 

z(t) = --W(e pt ) = -t+-\nW(e pt ) = -\n\e- pt W(e pt )] , (4.8) 
P P P 

where the equivalence between the expressions follows from the property lnW^z) = lnz — W(z). We refer to this 
trajectory as Proex, which is short for Productlog Exponential. A plot of its behavior is overlaid in Fig. [T] The late 
time behavior of this mirror is similar to the early time behavior of the Carlitz-Willey trajectory, in that it approaches 
the speed of light at timelikc infinity while not producing (in this case) a future horizon. The behavior is best seen in 
the Penrose diagram in the right panel of Figjl] Mathematically it can be seen by noting that as x — > 00, W(x) 00. 
Consequently, the value of v for the mirror at a given time, v m (t), has the behavior v m (t) — > 00 as t — > 00. The 
velocity is 

z{t) = [1 + W{e pt )] _1 - 1 , (4.9) 
which makes obvious the approach to lightspeed as t — > 00. The proper acceleration, 

<* * = -P [ ' , /2 , (4-10) 

[1 + 2W(e pt )f 2 

is initially zero, increases with time until it reaches a maximum magnitude when the trajectory intersects the null ray 
v = 0, then decreases with time, vanishing in the limit t — > 00. 

For this trajectory both the ray-tracing function p(u) and its inverse f(v) can be computed analytically. The results 
are 

p(u) =u- -W(2e KU ) , (4.11a) 

K 

f(v)= 2 e KV +v. (4.11b) 

K 

The Bogolubov coefficients /3 UW < can in turn be calculated analytically, as well as the energy flux. We summarize 
these and some other quantities in Table |III| 



D. Modified Carlitz-Willey trajectory 



The final class of trajectories that we consider is a modification of the Carlitz-Willey trajectory. A term is added 
that takes the acceleration to zero at late times. This causes the trajectory to become inertial in the future rather 
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TABLE III. Bogolubov coefficients and other analytically-derived information for the Proex trajectory. 



i(w+u/)/p 



Bogolubov coefficient 

Proper acceleration 

Spectrum of produced particles j^] 

Energy flux 

Total energy 



2-7ruip(ui+uj') e 2w(w+u')/P—l 



1 - e 



-27T(jj/ p 



alt) - -o W{ePt) 

,rp \ _ p 2 [2-W{2eP M )]W(2eP M ) 
\J-uu) — 4g7r [l + W(2eP«)] 4 



^E^!r(0,2^m/p) 



2-Kp 



Here r(0, 2iru>m/ p) is an upper incomplete gamma function. 



TABLE IV. Bogolubov coefficients and other useful information for the modified Carlitz-Willey class of trajectories. 



[Vj T[i{u + vu/)/k] 



Bogolubov coefficient 

Proper acceleration 
Energy flux 



2 



^7 e -7r(w+<rw')/(2 K ) 



a ® _ 2[<T+(l + t T)W( e - 2 «'/(l + -)/(l + t T))] 3/2 

Irp \ K 2 l-2fre K " 

\-Luu) — 487r ( CTeKU + 1 )2 ■ 



2'!TKijJ / (id + auj' ) e 27r(uj + (TUj / )/K_^ 
K ,(l + < 7)W( y £- 2Kt / {1 + "^ /(l + O-)) 



than asymptotically null. The final velocity is a free parameter. The trajectories are 

I- a 1 / e -2«t/(i+^) 
z = t- - W — 

K \ l + O 



1 



(4.12) 



with Q < a < 1. A particular example is shown in Fig. [T] where one can see the divergence from the Carlitz-Willey 
mirror at late times. For a — the trajectory reduces to the Carlitz-Willey one (3.4). As before with the other 
new trajectories, it proves possible to calculate the Bogolubov coefficients (3 uu i analytically, as well as a set of other 



relevant quantities. These are displayed in Table IV 



V. ANALYSIS OF ENERGY AND PARTICLE PRODUCTION 



For each of the mirror trajectories we consider in this paper, the renormalized stress-energy for the scalar field 
is known analytically and can be easily evaluated. Except for the Carlitz-Willey trajectory, the expectation value 
of the number of particles in a wave packet, (Nj n ), cannot be calculated analytically and we instead evaluate this 
quantity numerically. To do so accurately, however, it proved important to have analytic expressions for the Bogulobov 
coefficients j3 uu i . The wave packets in turn provide a means of examining the time and frequency dependence of the 
created particles. 



A. Time dependence 



The correlation (or lack of it) between the number of particles produced at a given time and the energy flux 
(T uu ) differs markedly from one type of trajectory to the other. Because the mirror is in flat space the Bogolubov 
transformation between the in and out vacuum states tells us about the average number of particles produced in an 
ensemble of identical systems. The energy flux given by the quantity (T uu ) gives information about the average flux 
of energy produced by the mirror as it accelerates. The energy flux is due to a combination of particle production 
and vacuum polarization effects. 

In Fig. [2] the energy flux is shown for the various trajectories considered in this paper. From Eq. ( |2.30[ ) one can 
see that the sign of the flux is closely tied with the change in the proper acceleration of the mirror. In particular 



14 




FIG. 2. The energy flux (T uu ) versus time is plotted for the various mirror trajectories. The parameters k, /j,, p, and p have 
all been set equal to 1. For the Darcx trajectory the value £ = 1/2 was chosen and for the modified Carlitz-Willey trajectory 
cr = l/3 was chosen. For the Walker-Davies trajectory A = 2 and B — 1 were chosen. The energy flux in the Carlitz-Willey 
case is the constant solid line. The energy flux associated with the modified Carlitz-Willey trajectory is the solid curve, which 
coincides with the Carlitz-Willey value at early times but then diverges, briefly resulting in a burst of negative energy before 
decaying to zero. The flux associated with the Arctx trajectory is shown as the short dashed curve and that of the Darcx 
trajectory is indicated by the long dashed curve. The flux from the Proex trajectory is depicted by the dotted curve and that 
of the Walker-Davies case by the dot-dashed curve. 



as mentioned in Sec. [TTJthe flux is negative if the change in the proper acceleration is towards the right and positive 
otherwise. 

Figs. [3] and [4] show the expectation value of the number of particles produced (Nj n ) as a function of the time 
parameter n for the various trajectories we consider. In Fig. [3j the frequency parameter j is set to 1 because 
of infrared divergences that occur in N 0n for the Carlitz-Willey and modified Carlitz-Willey trajectories. For the 
trajectories in Fig. [4j j is set to zero because no such divergences occur. The packets of course sample the particle 
production discretely. We draw attention also to the small level of particle excitation that occurs for j = 1 in Fig. [3] 
as compared to that which occurs for j = in Fig. |4j We return to this issue in Sec. |V C| 

There is a correlation between the number of particles created during a given time period and the flux of energy 
which occurs at that time for the Carlitz-Willey trajectory in Fig. [3] due to the fact that, as discussed in Sec. Ill A 



the flux is constant in time and the number of particles (Nj n ) is independent of the value of the time parameter n. 
As can be seen in Figs. [2] and 3 a correlation also occurs for the modified Carlitz-Willey trajectory at early times. 
However, the direct correlation is destroyed by the existence of a negative flux of energy after the time t = when 
the trajectory has deviated significantly from the original Carlitz-Willey trajectory. 

For the trajectories in Fig. [4] there is no direct correlation between the energy flux and the number of particles 
created. In fact for the Arctx trajectory there is something of an anti-correlation in that at about the time of peak 
particle production the flux is negative and has its greatest magnitude. This shows clearly the limitations in using the 
stress-energy tensor to describe the number of particles created. Because of vacuum polarization effects, which can 
include fluxes of negative energy, it is virtually impossible to separate out the contribution from the created particles. 



15 



3.0e-14 



2.0e-14 



1.0e-14 



O.Oe+00 



w ° o o o I 
_i i i i I i i i i I i i i i I i i i i 



-10 



-5 



10 



n 



FIG. 3. The particle number as measured by wave packets, (Nj n ), is plotted as a function of the packet time parameter n 
for the Carlitz-Willey (filled circles) and modified Carlitz-Willey (open circles) trajectories. In both cases the packet frequency 
width parameter e has been set to y/2-K and the first non-divergent frequency bin, j = 1, is shown. The parameter k has been 
set to 1 and for the modified Carlitz-Willey trajectory we have taken a — 1/3. 



B. Frequency spectrum and simultaneous frequency and time resolution 



Wave packets can also be used to measure the frequency spectrum of the created particles. Good frequency 
resolution is obtained by choosing a small value for e. As an example consider Fig. [5] for the = 1 Arctx trajectory 
with a wave packet frequency width parameter e = 0.01. The spectrum is a function of packet index j, and we have 
set n — 0. Clearly good frequency resolution is obtained. 

In principle the wave packet formalism allows particle production to be simultaneously resolved, within limits, in 
both frequency and time. The limits, of course, are set by the uncertainty relation. For any chosen e, the wave packets 
have a width in frequency Auj = e and an effective gating interval (width of time over which the particle detector 
is on) of At ~ 2tt/c. Thus, the time-bandwidth product (for these packets) is 2ir. It may be possible to resolve 
time-dependent spectra for a process with sufficiently copious particle creation and for which changes in its spectrum 
occur over a long enough time interval. However, for the trajectories studied in this paper, it has not been possible 
to resolve the creation process in frequency and time simultaneously, such that a significant number of particles is 
found in many frequency and time bins. Some of the mirror trajectories have confined periods during which the 
acceleration peaks and is strong (i.e., Walker-Davies, Arctx, Darcx, and Proex). For these trajectories, when we 
choose a relatively small value of e to provide good frequency resolution, we find the vast majority of the particles are 
created in the time bin n = 0. In contrast, when we choose a relatively large value of e to gain good time resolution, 
we find that almost all of the particles reside in the lowest frequency bin, j = 0. Furthermore it is not even possible 
to find some intermediate value of e for which some time and frequency resolution is possible. Instead, in picking such 
an intermediate value of e, we find that the vast majority of particles lie in the single bin with n = j = 0. 

A different behavior occurs for the Carlitz-Willey trajectory. In this case a small e will provide adequate frequency 
resolution and yet not all of the particle creation occurs in a single time bin, such as n = 0. However, we have not 
really succeeded in simultaneous time and frequency resolution, since the Carlitz-Willey mirror gives a spectrum that 
is completely time independent. This then brings us to the modified Carlitz-Willey trajectory. The energy flux in 
this case is asymptotically constant in the distant past but then at some point the energy flux drops toward zero 
as the acceleration falls off and the mirror becomes inertial. The particle creation behaves similarly. One might 
hope that this trajectory would result in a creation process that could be simultaneously resolved in frequency and 
time. However, here too we find that the acceleration falls off sufficiently rapidly that the transition from creation to 
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FIG. 4. The particle number as measured by wave packets, (Nj n ), is plotted as a function of the packet time parameter n. 
In each case the packet frequency width parameter e has been set to y/2ir and the lowest frequency bin j = is shown. The 
parameters fi, v, and p have all been set to 1. For the Darcx trajectory the value £ = 1/2 was chosen and for the Walker-Davies 
case A = 2 and B = 1 were chosen. The value of No n is denoted by the open circles for the Arctx trajectory, the triangles for 
the Darcx trajectory, the filled circles for the Proex trajectory, and the pluses for the Walker-Davies trajectory. 



effectively no creation occurs within one time bin (assuming e has been set to allow good frequency resolution, i.e. 
many frequency bins within the transition or characteristic frequency u> c ~ k) . 

To see what seems to be happening, consider the Arctx trajectory. In that case the parameter fi is dimensionally 
the inverse of time and ought to represent a characteristic frequency u c . In fact, as Fig. [5] shows, we find the peak 
of the particle creation spectrum is uj c ~ 0.14/1. However, the period r of significant acceleration also depends on 
[i and is roughly r ~ \i . The frequency and time scales are thus related by a single parameter. A related factor 
is that the expectation value of the total number of particles created (N) is a fraction of unity. In other words, in 
an ensemble of identically-accelerated mirrors, in many cases there will be no particles produced at all. Similar time 
scale and frequency scale issues occur with the Darcx, Proex, and Walker-Davies trajectories. 

We speculate that a trajectory might be crafted with two parameters: an acceleration scale /z, with r ~ [i _1 , and 
a duration of acceleration T, which satisfies T 3> r. For such a trajectory enough particles may be created for a 
long enough period to allow a significant number of particles to be found in many time and frequency bins, giving 
good resolution. A mirror of this sort would undergo a large change in Lorentz factor over a time ~ T and it would 
appear, over that time interval, like a mirror that is approaching a null horizon. We have not been able so far to find 
a trajectory with these properties for which the Bogolubov coefficients can be calculated analytically. 

The challenge of finding significant wave packet excitation that is simultaneously spread across a range of both 
time and frequency bins may also lie in the inherent nature of the time-bandwidth product of our orthonormal wave 
packets. The time-bandwidth product of these packets is 27r, while the fundamental limit of the uncertainty principle 
is 1/2. This may be a contribution to non-uniformity in the creation spectra, since by their construction, these 
orthonormal wave packets are unable to reach the limits of the uncertainty principle. 

The fact that it is not possible to obtain significant particle creation in simultaneous bins of both frequency and 
time for the trajectories considered may have interesting experimental consequences if a system that was in some 
way like one of these trajectories could be studied in the laboratory. However, the relatively small amount of particle 
production that occurs might make this very difficult. 
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FIG. 5. The particle number as measured by wave packets, {Nj n }, is plotted as a function of the packet frequency parameter 
j, with n = and e = 0.01 for the Arctx trajectory with (M= 1. 



C. Total number of particles produced and their energy 



In Sec. [TT| expressions were given that allow one to compute the total number of particles produced and their total 
energy. With some mirror trajectories, we have found one or both of these quantities to be divergent. 

To compute the total number of particles produced one can use Eq. (2.221 and to compute their total energy one 
can use either Eq. (2.261 or (2.31). Comparison of Eqs. (2.21), ( |2.22 1, and (2.26) shows that it is possible to have a 



divergent number of particles produced and yet have a finite total energy, provided that the divergence in the number 
of particles is due to an infrared divergence in N u and that it is not too strong. 

Substitution of |/3 WW '| 2 into (2.22) and (2.26) for the Walker-Davies and Arctx trajectories results in both a finite 
number of particles produced and a finite total energy for those particles. We have evaluated (N) numerically for 
both trajectories. For the Walker-Davies trajectory, (N) depends on the ratio < B/A < 1. As an example, for 
B/A = 1/2 we find (N) rj 0.0121. For the Arctx trajectory, (N) is independent of the value of the parameter /i 
as can easily be seen by substituting (4.5) into (2.22) and making variable changes of the form u> — fix. We find 
(N) « 0.0134. An interesting point is that both of these values are small compared to unity. Thus in an ensemble of 
identical systems, in most cases there would be no excitation of any mode and no particles would be produced. The 
total energy produced for the Walker-Davies trajectory was computed by them and is given in Eq. (3.13). For the 
Arctx trajectory it is given in Eq. (4.4). As discussed in Sec. [ill the energies E st in (2.31) and E qs in (2.26) are the 
same. This can be used as a check on the computations of the Bogolubov coefficients P^'. We have computed E qs 
numerically (in specific cases for the Walker-Davies trajectories and in general for Arctx) and found the values to be 
equal to those for E st . 

For the other classes of trajectories besides Arctx and Walker-Davies, the number of particles (N) diverges. For 
the Proex trajectory we have computed (JV ( 
indicates an infrared divergence of the form 



analytically. Examination of the result, which is shown in Table III 



ln(27rw/p) 



UJ 



(5.1) 



This behavior will result in a divergence in the total number of particles (N) . However when computing E qs in 



Eq. (2.26) one multiplies {NJ) by a factor of co before integrating over uj. The result is a finite value for E qs . We have 
computed E qs numerically and found agreement with the expression for E st , which is given in Table III 
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For the Darcx class of trajectories when u is small, one can divide the integral in ( |2.21 ) into two parts, one with an 
integral from to v\ and a second integral from v\ to oo with 1 > A > u/v. In the first integral one can evaluate 
|/3 wa ,'| 2 in the limit that both u/v and u' /v are small. The result is a contribution to (N u ) that is proportional to 
1/uj. In the second integral, u <C u' and it is not hard to show that the integral is finite in the limit u — > 0. Thus, as 
with the Proex class of trajectories, we find that (N) is infinite but E qs is finite. We have numerically computed E qs 
for specific values of the parameter £ and shown that its value is the same as that for E st , which is given in Table |H| 

For the Carlitz-Willey trajectory, one can see from the form of |/3 WW '| 2 that there is both an infrared and ultraviolet 
divergence in the integral over u' in Eq. (2.21 ). Thus there is a divergence in both the particle number and the energy 



of the produced particles. The latter is trivially apparent from the constant flux of energy which occurs f or t hese 

v I 2 in Table: 



IV 



that 



0. However, for all other packets (Nj n ) is finite so the divergences are 



trajectories. For the modified Carlitz-Willey class of trajectories one can see from the form of \/3 u 
there is still the infrared divergence but no ultraviolet divergence when computing (N u ). The infrared divergence is 
strong enough to make (N u ) divergent for all values of u. 

One might ask what effect the packets have on those trajectories with infrared divergences. For the Carlitz-Willey 
trajectory we were able to analytically compute the packets and as seen in Eq. (3.8) there is a divergence for the 
packets with the lowest frequency range, j 
not nearly as strong as for (iV w ). 

For the modified Carlitz-Willey class of trajectories there is also a divergence for packets with j 
those with j > 0. To see this, one can divide the integral in (2.38) into two parts such that Ii = f*" du' and 
I2 = / » du' , with < A <C 1. For the second integral, which contains only nonzero values of u', it is not difficult 
to see that for f3j nM > in Eq. (2.37) there are no infrared divergences resulting from the integral over u. It can also 
be seen that fijn^ 1 is we U enough behaved in the limit u' — > 00 that there are no ultraviolet divergences so long as 
a > 0. For the first integral, the value of u' / k is small. Here it is necessary to break the discussion into the cases 
j = and j > 0. For j > 0, we take a small enough value of A so that e > u' for all values of u' in the first integral. 
Then one can expand the function r[z(u; + u')/k] in the expression for /3 wt y in Table IV in powers of uj' /uj. Keeping 
the leading order term one finds that 



but not for 

A 



\0. 




^V^e [2wm/£ - 7r/(2K)I " e -i<W«)M"7*) T{iuj/K) 



(5.2) 



By repeatedly integrating by parts one can obtain a series in inverse powers of ln(u/ /k). Substituting into 1\ then 
shows that to leading order the integrand goes like l/[u/(ln(u///i)) 2 ], which when integrated gives no divergence in 
the limit u' — > 0. 



For j = the situation is different. Here one can divide the integral in (2.37) into two parts, J\ — L duo and 
J2 = f » duo. The analysis for the second integral is exactly the same as for the case j > 0. For the first integral we 
choose A to be small enough so that \2irKn/e — 7r/2|A <C 1. Then both u/k and u' / k are small and 




duJue- l(uj/K)ln(ul ' /K) - 



u + u' 



(5.3) 



Making the change of variable z = —(u/k) \n(u' /k), the upper limit becomes —XIh(u'/k), which for fixed A goes to 
infinity in the limit u' — > 0. The resulting integral can be computed analytically in terms of Fresnel Integrals. The 
result to leading order in u' is Ji - l/y/-u'ln(u'/K). This gives a contribution to I\ that when integrated results in 
a divergence at the lower limit u' — and thus a divergence in (Nj n ) for j = 0. 

For the Darcx and Proex trajectories it turns out there is no infrared divergence in (Nj n ) even for j ' = and even 
though (NJ) diverges in the limit u — > 0. This can be shown by bounding (Nj n ) by substituting Eq. (2.37) into 
Eq. (2.38) and then computing the absolute values of each factor in the integrand of (2.37). The result is 




(5.4) 



It is not hard to show that for small u' the integrand is finite in the limit u' — > if j > for both trajectories. If 
j = then it is useful to again divide the integral over u into two parts as was done above, except with k — ¥ v and p 
respectively for the Darcx and Proex trajectories. Then as before the analysis for Ji is the same as for j > 0. For Ji 
both u and u' are small making it possible to expand the terms in |/3 WW '|. When this is done and the integral over u 
is computed for the leading order terms, we find that the result is finite in the limit u' — > for the Darcx trajectories. 
For the Proex trajectory, we find that to leading order Jj ~ In a/ so that the integrand for the integral over u' goes 
like (In a/) 2 and a finite contribution is made to (No n ). 

As discussed in Sec.[TTJ one can obtain an estimate of the energy of the created particles using the packet formalism 
by multiplying (Nj n ) by the frequency in the middle of the range for each packet and summing over j and n as in 
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Eq. (2.39|. The resulting energy, E ep , has been computed for the Arctx trajectory for two different values of e. For 
the case shown in Fig. [5] with e = 0.01 the results agreed with E st in Eq. (4.4 1 to within about 0.01%. In a separate 
calculation, with e = 10, the agreement was at the 1% level, which is remarkably good given the poor frequency 
resolution, which might be expected to drastically skew the energy summation. 

For the Darcx class of trajectories with £ = 0.99, a = 1, and e = 0.001, we find agreement with the values of E st 
to within four digits. The result was obtained by summing packets with n = and values of j ranging from j = to 
3 = 1175. 

For the Proex trajectory with p — 1, a numerical computation of E ep for e = 5x 10~ 5 gave results in agreement 
with E s t to approximately 0.08%. Energy packets with n — and with a sum from j = to j = 12000 were used in 
that calculation. 



VI. CONCLUSION 



We have investigated the particle production and the energy flux that results from a massless, minimally coupled 
scalar field in a two-dimensional flat space that contains an accelerating mirror. Dirichlet boundary conditions are 
assumed at the mirror and the field is assumed to be in the in vacuum state. Six different types of trajectories have 
been considered, including the one studied previously by Carlitz and Willey [7] and the one studied by Walker and 
Davies [22]. The other four are new and have been introduced for this study. These trajectories are all asymptotically 
inertial in the limit t —00 and all but the Carlitz- Willey trajectory are also asymptotically inertial in the limit 
t — > 00. 

For each trajectory it has been possible to obtain analytically the Bogolubov coefficients /3 UU > as well as the proper 
acceleration and the energy flux (T uu ). As pointed out by Walker [IT], it is very useful to have models in which it is 
possible to do analytic calculations. The four new types of trajectories that we have provided fit this description. 

Our main focus has been on the use of wave packets, which allow the particle production to be time resolved, and 
in principle might allow significant simultaneous frequency resolution as well. The packets we use form a complete 
orthonormal set, so that no information is lost. By computing the Bogolubov coefficients for the packets and inte- 
grating over the frequencies of the in modes, it is possible to obtain the average number of particles reaching 
found in a given frequency range and an approximate time range (more specifically a range in the null coordinate u) . 
Thus this method of analysis can be thought of as similar to what a series of particle detectors along a large v surface 
would detect if each was turned on for some relatively short period of time. 

In principle one might expect a correlation between the time dependence of the particle production and the energy 
flux (T uu ). However both vacuum polarization and particle production effects are combined in the stress-energy tensor 
and it is difficult if not impossible to separate them. Use of the wave packet formalism allows for an unambiguous 
description of the time dependence of the particle production process. 

The Carlitz- Willey trajectory was designed to result in a constant flux of energy and is of a different nature than 
the other trajectories in being asymptotically null. We found an explicit mathematical expression for this trajectory 
in terms of the Lambert W function. Because of the constant flux of energy the total energy produced is divergent. 
Not surprisingly it is also found that the number of particles produced per frequency interval (N^) is also divergent. 

This trajectory was the only one for which we were able to compute the number (Nj n ) of particles associated with 
a wave packet analytically. Just as the energy flux is constant, we found that (Nj n ) is independent of the value of n 
and so is time independent. A divergence occurs for the case j — but not for larger values of j. This infinity can be 
dealt with by simply ignoring the lowest frequency (j = 0) bin. In a realistic detector there will always be an infrared 
cutoff, since it is impossible to detect particles of arbitrarily long wavelengths. An exploration of the effects on the 
frequency range for the wave packets with j > was carried out and it was found that a Planck type spectrum is 
approached in the limit that the frequency width of the packets vanishes. 

The modified Carlitz- Willey trajectory has the same approximate behavior as the Carlitz- Willey trajectory at 
early times but then stops accelerating and approaches a constant velocity at late times. As with the Carlitz- Willey 
trajectory the number of particles produced per frequency interval (iV w ) is divergent for all ui. The energy flux (T uu ) 
is approximately constant at early times and this results in an infinite amount of total energy E st . Using wave packets 
we again find that (Nj n ) is divergent for j = but finite for all other values of j. For j > we find that, as a function 
of the time parameter n, (Nj n ) is approximately constant at early times and decreases to zero at late times as would 
be expected for a trajectory that is asymptotically inertial. 

The Arctx and Walker-Davies trajectories are the only ones for which the mirror begins and ends at rest. We find a 
finite number of particles (N) is produced for both along with a finite amount of energy E st . Using wave packets we 
have shown that the number of particles produced (Nj n ) increases to a maximum and then decreases over the range 
of time that the mirror's acceleration is first increasing and then decreasing at a significant rate. 
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For the Darcx and Proex trajectories the mirror begins at rest and is asymptotically inertial in the future. In the 
Darcx case it approaches a constant speed that is less than that of light and in the Proex case it approaches the speed 
of light but in such a way that it is not asymptotic to a null trajectory. For these trajectories there is an infrared 
divergence in the number of particles produced but not in the energy of the produced particles. However there is no 
divergence in (Nj n ) for j = 0. We find that the number of particles first increases to a maximum and then decreases 
during the period when the acceleration is first increasing and then decreasing at a significant rate. 

It is interesting to compare the results for the time-dependent particle production with the average energy flux (T uu ). 
Not surprisingly, for the Carlitz-Willey trajectory both are constant in time, so in that sense there is a correlation. 
For the modified Carlitz-Willey class of trajectories there is a similar correlation at early times. Once the number of 
particles produced begins decreasing, the correlation diminishes and even disappears due to a flux of negative energy 
that occurs at intermediate times. For the other trajectories, probably again because of the fluxes of negative energy 
during certain time periods, there is little or no correlation between the number of particles produced and the energy 
flux. 

Because the wave packets tile both the time and frequency domains, it might be expected that one could obtain 
time-dependent spectra for the particle production. For each type of trajectory (except Carlitz-Willey) the time 
and frequency ranges of the bins were varied. It was found (at least for these trajectories) that simultaneous time 
and frequency resolution with significant particle build-up in each domain is absent. The mirrors have a single 
dimensional parameter that determines both the characteristic frequency and duration of creation, and for this reason 
the uncertainty principle prevents one from measuring the spectral dynamics. It may be possible to find mirror 
trajectories with two characteristic scales, one that sets the acceleration and characteristic frequency and one that 
sets a duration of creation. This would allow significant particle creation in the spectral dynamics to be measured by 
wave packets. We have not found such a trajectory. Conversely, it may be a generic feature of the quantum nature 
of accelerating mirrors. Our results are in two dimensions, not four. However, it seems unlikely that this effect of the 
uncertainty principle is tied to the number of dimensions. If our results are pointing to a generic effect, then it could 
have important observational consequences for any experiments that might attempt to detect the radiation produced 
when a mirror accelerates. 
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